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Unser DUG-Mitglied Johannes Zerbs ersucht um Bekanntgabe seiner Mathe Chilled 
Website. Dieser Bitte komme ich gerne nach. Was besonders bemerkenswert ist: Es 
handelt sich hier um kein kommerzielles Angebot! Ja, so was gibt es auch noch, Josef. 



Mathe Chilled- ein Online Mathe Buch: 

Ein Ratselblock nach den neuesten Standards und Kompetenzen 

www.mathe-chilled.at 


Gestem Mathe. Heute Mathe Chilled. Morgen in Mathe chillen.... 


As TIME Goes by ... 


Past TIME 


and future TIME 




The “Kissing Students” are one of the most famous landmarks of Tartu, hosting town of 
TIME 20 12 (left). 

View of the old town of Krems (right). Krems is the entrance to the Wachau which is 
surely the most romantic part of the Danube valley between the Blackwood Forest and 
the Black Sea. Krems will host TIME 2014. Here is the place where it all began in 1992 
with the first DERIVE Conference. 
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Liebe DUG-Mitglieder, 

leider bin ich wieder einmal nicht 
rechtzeitig mit unserem Newsletter fertig 
geworden. 

Einerseits war ich sehr mit anderen 
Aufgaben beschtiftigt und andererseits war 
dieser DNL sehr zeitaufwandig. Die Uber- 
setzung des Statistikartikels ging rasch 
von statten - ermdglicht durch die gute 
Kommunikation mit Guido Herweyers. Dage- 
gen brauchte mein eigener Branch <5t Bound- 
Beitrag viel Zeit bis er meinen Wunschen 
entsprach. Aber das war ein schoner Aus- 
flug in meine Lehrervergangenheit. 

Piotr Trebisz hat den letzten Teil sei- 
ner Schneckenhausserie geschickt. Ein 
spannender Austausch von E-Mails ermo'g- 
lichte interessante Erganzungen zur Dar- 
stellung seiner schonen konischen Spiral- 
flachen. Aus Platzgriinden muss ich das fur 
den ntichsten DNL zuruckstellen. 

Auch Dietmar Oertel wartet schon lan- 
ge auf die Veroffentlichung seiner Arbei- 
ten an den Primzahlen und Taylorreihen. 
Ich ersuche auch ihn um Geduld. Sie sind 
„under construction". 

Noch kurz zur TIME 2012: Herzlichen 
Dank an die Organisatoren vor Ort, Marina 
und Eno samt ihrem Team an der Universi- 
ttit Tartu, Estland. GroGen Dank auch an alle 
Teilnehmer, die hervorragende Hauptvor- 
trage und hochklassige Lektionen angebo- 
ten haben. 

Eine Reihe von Vortragenden hatte 
nicht die Absicht, ihre Beitrage zur Verof- 
fentlichung in Journalen einzureichen. Sie 
haben ihre Papiere dem DNL zur Verf iigung 
gestellt. Vielen Dank dafur. 

SchlieBlich mdchte ich noch besonders 
auf die Themen im User Forum aufmerksam 
machen. Sie finden zwei kurze Hinweise auf 
TT-Nspire-Funktionen und zwei Probleme, 
die auf unsere DUG-Mitglieder warten. 

Ganz aktuell: gemeinsam mit der Donau 
Universittit Krems wird ACDCA die TIME 
2014 vom 1. - 5. Juli 2014 in Krems an der 
Donau veranstalten. 


Dear DUG Members, 

unfortunately I am too late with this 
newsletter once more. 

On the one hand I was busy with other 
duties and on the other hand production of 
this DNL was very time consuming. Transla- 
tion of the statistics contribution was - 
thanks good communication with Guido 
Herweyers - quickly done. My own article on 
Branch A Bound took much time until it met 
my expectations - and during writing new 
ideas wanted to be realized ... This was a 
nice trip back to my past as teacher. 

Piotr Trebisz sent the last part of his 
snail shell series. An exciting email ex- 
change brought interesting additions to the 
representation of his beautiful conic spiral 
surfaces. Because lack of space I must 
leave this for the next issue. 

Dietmar Oertel has also been waiting a 
long time for publication of his prime num- 
ber and Taylor series investigations. Please 
be patient. They are "under construction". 

Concerning TIME 2012: many thanks to 
the local organizers, Marina, Eno and their 
team at Tartu University, Estonia. Many 
thanks also to all participants and present- 
ers. We had excellent keynotes and great 
lectures. 

Some presenters didn't intend to sub- 
mit their full papers for publication in spe- 
cial journals. They provided their presenta- 
tions for the DNL. Thanks for this offer. 

Finally I'd like to draw your attention to 
the User Forum. You will find two short 
notes on Nspire functionality and two prob- 
lems waiting to be tackled by DUG mem- 
bers. 

Latest - great - news: The Danube Uni- 
versity Krems and ACDCA will organize 
TIME 2014 from 1 to 5 July 2014 in Krems, 
Lower Austria. 

Viele GruBe, kindest regards 



Download all DNL-Derive- and TI-f!les from 

http : / / www . austromath . at/ dug/ 
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The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE & CAS- 7/ User Group. 
It is published at least four times a year 
with a content of 40 pages minimum. The 
goals of the DNL are to enable the ex- 
change of experiences made with DERIVE, 
TI- CAS and other CAS as well to create a 
group to discuss the possibilities of new 
methodical and didactical manners in 
teaching mathematics. 


Editor: Mag. Josef Bohm 
D'Lust 1, A-3042 Wiirmla, Austria 
Phone: ++43-(0)660 3136365 

e-mail: nojo.boehm@pgv.at 


I II 

Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
the DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE & CAS- 77 Newsletter will be. 

Next issue: December 2012 


Preview: Contributions waiting to be published 

Some simulations of Random Experiments, J. Bohm, AUT, Lorenz Kopp, GER 
Wonderful World of Pedal Curves, J. Bohm, AUT 
Tools for 3D-Problems, P. Liike-Rosendahl, GER 
Hill-Encryption, J. Bohm, AUT 

Simulating a Graphing Calculator in DERIVE, J. Bohm, AUT 

Do you know this? Cabri & CAS on PC and Handheld, W. Wegscheider, AUT 

An Interesting Problem with a Triangle, Steiner Point, P. Luke-Rosendahl, GER 

Graphics World, Currency Change, P. Charland, CAN 

Cubics, Quartics - Interesting features, T. Koller & J. Bohm, AUT 

Logos of Companies as an Inspiration for Math Teaching 

Exciting Surfaces in the FAZ / Pierre Charland 's Graphics Gallery 

BooleanPlots.mth, P. Schofield, UK 

Old traditional examples for a CAS - what’s new? J. Bohm, AUT 

Truth Tables on the TI, M. R. Phillips, USA 

Where oh Where is It? (GPS with CAS), C. & P. Leinbach, USA 

Embroidery Patterns, H. Ludwig, GER 

Mandelbrot and Newton with DERIVE, Roman Hasek, CZK 

Tutorials for the NSpireCAS, G. Herweyers, BEL 

Some Projects with Students, R. Schroder, GER 

Structures in the Set of Prime Numbers, D. Oertel, GER 

Dirac Algebra, Clifford Algebra, D. R. Lunsford, USA 

Laplace Transforms, ODEs and CAS, G. Picard & Ch. Trottier, CAN 

A New Approach to Taylor Series, D. Oertel, GER 

Cesar Multiplication, G. Schodl, AUT 

Henon & Co; Find your very own Strange Attractor, J. Bohm, AUT 

Rational Hooks, J. Lechner, AUT 

Mathematical Model for Snail Shells (4), P. Trebisz, GER 

Simulation of Dynamic Systems with various Tools, J. Bohm, AUT 

An APL-like SHAPE function in DERIVE 6, D. R. Lunsford, USA 

Brussels Gate in Dendermonde, E. van Lantschoot, GER 

and others 


Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wiirmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 
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Statistics with Tl-Nspire 3.1/3.2 (Part 3) 

Visualising and Simulating Dynamically with Tl-Nspire 

Guido Herweyers, KHBO Campus Oostende 
quido.herwevers@khbo.be 


Part 2: Simulating with Tl-Nspire 

(1) Generating random numbers 


Open a calculator page. In the menu 5: Prob- 
ability > 4: Random can you find the com- 
mands for generating random numbers. 

Click on Seed or type randseed followed by 
space. Then type any natural number in order 
to initiate a new series of random numbers. 


Examples for generating random numbers 
(you may compare with the TI-84 Plus com- 
mands) are: 


Documents Toolbox 


f>i 


Actions 


X = 3: Algebra 
Id 4: Calculus 


5: Probability 


J 6: Statistics 

7: Matrix Si Vector 


010 


^Functions Si Programs 


1 : Factorial (!) 

2:Permutations 

3:Combinations 


5: Distributions 


1 :Number 
2: 1 nteger 
3:Binomial 
4:Normal 
5:Sample 
6 : S e e d 


rand(): 

rand(IO): 

randint(1,6): 

rand(1 ,6,20): 

randbin(10,0.5): 

randbin(10,0.5,20): 

randnorm(175,10): 


a real random number between 0 and 1 . 

So gives 3 + 5 • rand() a random number within the interval (3,8). 

a list of 10 random numbers between 0 and 1 . 

a natural number from 1 to 6 (throwing a die). 

a list of 20 integers within [1,6] (throwing a die 20 times). 

the number of heads when throwing a coin 10 times. 

a list of the results repeating “= TIMES the coin experiment from above. 

a normal distributed random number with mean =175 and 
standard deviation =10. 


randnorm(1 75, 1 0,50): a sample (list) of 50 random numbers from a normal distribution with 

mean =175 and standard deviation =10. 


RandSeed 834 

Done 


randlnt(l,6,10) 

{ 5, 5, 6, 3, 5, 5, 5,5, 1,5 } 


randlnt(l,6,10) 

{2,3,1,3,4,2,6,6,1,31 


randBin( 10, 0.5,20) 

{ 3, 2, 7, 6, 6, 3, 4, 7, 4, 6, 6, 6, 6, 5, 6, 3, 7, 6, 4, 6 } 


3+5- rand (4) 

{ 4. 84643,6. 80426,7. 00736,4. 05879 } 
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(2) Samples with and without replacement 

A sample from a list of numerical or categorical data is taken using the command randsamp. Its syn- 
tax is: 

with replacement: randsamp(list name, sample size, 1) 

without replacement: randsamp(list name, sample size) 

Drawing a random sample without replacement: 
rand samp (Li sc name, sample size.1) 

Example: Lottery drawing (6 out of 45) 

RandSeed IDO * Done 
lottery =seq(^,^l J 45) 

rarLdSampflotter yAl) * {42,29,7,23,26,45} 

Drawing a random sample with replacement: 
rand 5 ami p (List name, sample size) 

Example: 10 turns of a roulette wheel 
r ovt ett e : =^seq £, 0 S 36) 

- { 0 r 1 A 3, 4, 5, 6, 7, S, 9, 1 0 , 1 1 , 1 2, 1 3, 1 4, 1 5, 1 6. 1 7 r 1 S, 1 9, 20,21, 22 r 23, 24, 25, 26, 27,28, 29, 30, 3 1, 32,33, 34, 35, 36 } 
randSaiTip(eDuittte, 10) - { 27, 32, 11,18,1 4, 5, 1 7, 24, 30, 2 1 } ] 

(You can find the roulette rules together with the win chances in 
http://www.casino-gids.be/artikels/spelregels/beginners/roulette.php.) 

The screen shown above was created in a Notes Application. 

The output of a Maths Box used in a Notes Application 
can be shown or hidden. 


Perform a right mouse click in the Maths Box and choose 
the respective attributes. 

The lottery list {1,2, . . ., 44, 45} is hidden and the 
roulette list is shown. 


Input & Output: 
Insert Symbol: 
Display Digits: 


Hide Output 




Show Input & Output 
Hide Input 


Hide Output 


Angle . |lMo Calculation 
0 Wrap expressions 
0 Show warning indicator 


The Notes Application is very useful. It is an elementary text processor combined with a dynamic 

calculator in an environment which enables easy quick repetitions of simulations. 

• You are able to combine text and mathematical expressions (like in a calculator page). A mathe- 
matical expression is written in a Maths Box, which can be activated by pressing Ctrl + M. As 
shown above you can hide the result of the evaluation: right click on the Maths Box and then 
choose the Maths Box Attributes. 

• The Notes Application is dynamic (in contrary to the static Calculator Application): changes of a 
definition result in an immediate recalculation of all expressions in the Notes. 

• A simulation - which is an expression containing rand... in any form - is automatically reevalu- 
ated by clicking on the expression followed by ENTER (or by pressing Ctrl+ENTER in order to 
perform quick repetitions of the simulation). Unfortunately this does not work any longer in ver- 
sion 3.2. Here you need a respective column in the spreadsheet where you can repeat the simula- 
tion by pressing Ctrl+R. Read more about this on page 6. 
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Hint: Having defined a column in a spreadsheet application using a simulation command then the 
simulation is immediately reevaluated after clicking the spreadsheet and then pressing Ctrl+R. 



Altimes 

Broils 

♦ 

| 

|=randint(1,6,a1 \ 

1 

10 

6 

2 


6 

3 


1 

4 


4 

5 


3 

6 


5 

7 


1 

8 


3 

9 


6 

10 


2 


Ctrl+R 



Altimes 

W rolls 

♦ 


=randint(1,6,a1) 

T 

10 

4 

2 


3 

3 

1 

1 6 

4 


6 

5 


4 

6 


5 

7 


4 

8 


5 

9 


2 

10 


2 


Ctrl+R 



^Itirnes 

drolls 

* 

T 


=randint(1 i 6 i a1) 

10 

4 

2 


5 

3 

1 

1 2 

4 


1 

5 


3 

6 


6 

7 


3 

8 


4 

9 


1 

10 


4 


(3) Rolling two Dice 

We split the screen into three applications: Lists & Spreadsheet, Notes, and Data & Statistics. Sample 
size n and simulation sumd (sum of the two dice) are defined within the Notes. Spreadsheet and bar 
chart as well will change dynamically when performing a new simulation by changing the sample size 
n and by repeating the simulation or by clicking in the sumd Maths Box and then pressing ENTER. 
Investigate the influence of the sample size on the variability of the sample and the “convergence” to 
the theoretic probabilities of the sum of two dice. 



A sumc j ISresu... 

C||f e q l®relfreq 

E probabilities 

F n 

* 

T 







9 

6 

12 

0.12 

0.1389 


6 

7 

7 

20 

0.2 

0.1667 


7 

4 

8 

13 

0.13 

0.1389 


8 

5 

9 

12 

0.12 

0.1111 


9 

9 

10 

6 

0.06 

0.0833 


10 

4 

11 

4 

0.04 

0.0556 


11 

7 

12 

4 

0.04 

0.0278 


12 

9 






13 

11 






14 

7 






15 

4 






16 

10 






17 

10 






18 

7 






19 

8 






20 

6 






21 

6 





a 

/ 5 \ 

E7 =round — f 4\ 
\ 36 j 







Simulation of rolling n times two dice and 
observing the sum. 

Rolls n: = 10G 

The sums of the dice is put in column A. 
sumd : =rand int( 1 , 6, n)+rand int( 1 , 6, n) 

The frequencies of the sums is put in column C. 
fr eq: =seq(countif(sumd J ? =k),k J 2, 12) 

The relative frequencies are given by 
relfreq:=approx(freq/n) . 


Compare these values with the resepctive 
probabilities! 


24% - 
1 8 %- 


0 




1 2 %- 
6 %- 
0% 


1 2 3 4 5 6 7 8 9 10 11 12 13 

sumd 


It is really a pity that the full text in this paper is only static and not dynamical . . . 
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(4) Probability distribution of the sample mean 

In the following we will discuss the principle how to find an approximation for the probability model 
of the sample mean. 

• Take a big number of samples from a population with a discrete or continuous probability 
distribution. 

• Note for each sample its mean and append its value to a list via automatic capturing in a 
column of the spreadsheet application. 

• A dot plot of the list shows whether the probability model of the sample means should 
be a discrete or a continuous one. 

• Observe graphically the convergence of the percental distribution of the bar chart for 
a discrete model or the convergence of the density histogram (with small enough class 
widths) for a continuous model to a stable situation. 

• The “equilibrium situation” gives an approximation of the probability distribution of the 
sample means. 

For didactical reasons it is recommended to take the sample from a concrete given population of num- 
bers like the results of main part 4 of the life guards which is a population of 367 data showing a nega- 
tively (= left) skewed distribution (long lower tail). 

This is much more realistic for the students than taking a sample from an artificial population which is 
generated by applying a probability model (e.g. the uniform distribution for rolling a die). 

Introduce the data as shown below in the Notes. Sometimes it is necessary to execute a command with 
Ctrl+ENTER (instead of only ENTER) in order to get the result in decimal form (instead of a frac- 
tion). Examples for this are samplemean and mean(h4). 

The mean of each sample will be automatically added to the list means, which has been defined in the 
spreadsheet application: click into the formula cell and call Data > Data Capture > Automatic in the 
Spreadsheet Toolbox, then change capture(var,l) to capture(samplemean,l). 

For taking random samples provide a Maths Box (Ctrl+M) and type sample:=randsamp(h4,n). 
Another Math Box is necessary to calculate the samplemeans: samplemean:=mean(sample). These 
results are collected in the spreadsheet column means. To get a really good impression you need 
many samples, so keep pressing Ctrl+ENTER in the sample-Maths Box for repeating the process 
many times. Accomplish the Notes according to the screenshot given below . . . but . . . 

... inspecting this screenshot you will find another column sample in the spreadsheet. Why this? Un- 
fortunately TI-Nspire changed its behaviour from Version 3.1 to 3.2 as mentioned on page 4 in such a 
way that pressing Ctrl+ENTER in the Maths Box does not repeat the simulation. You have to define 
the sample in a column. Then click into any cell of the spreadsheet and press Ctrl+R. 

Perform a sufficient large number of samples. 

Don’t forget to change the RandSeeds or ask the students to take different RandSeeds otherwise they 
could get identical results - and this is will not give a “random like” impression. 
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from h4 with sample size n:=4 
RandSeed 1010 ► Done 

sample " { 70. ,54., 83. ,89. } 

sampleiiieaii:=mean(sample) " 74. 

Number of calculated sample means: dim(means) * 2055 
Compare the population mean and the mean of the sample means: 

mean(h4) ► 65.0027 mean(means) > 64.9751 

Compare the standard deviation of the population and the standard deviation of the sample means: 

stDevPop(h4) 

-=r+ — - > 9.04752 

Jn 


stDevPop(li4) ► 18.095 


stDevSamp(means) * 8.81557 



■ sample 

■ means 

eg 

♦ 

=randsamp('h4,'n) 

=capture('samplemean,1 ) 


1 

70 

48. 


2 

54 

71.5 


3 

83 

67. 


4 

89 

54.25 


5 


65.25 


6 


64.75 


7 


| 66. 


8 


66.25 


9 


56.25 


10 


78.25 


11 


73.75 


12 


78. 


13 


54.25 


A 

"1 

r 

1 


0.08 


0.06 


0.04 


0.02 


0.00 



Working with sample size n = 4 we find the probability distribution of the mean less skewed than the 
populations’s one. You can see a simulation of 2055 samples. The average of all means is 64.98 which 
is very close to the mean of the population 65.00. The deviation is almost halved. 



Sample size n = 16 and a series of simulation of 2251 samples. Deviation is halved again. 
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from h4 with sample size n:=64 
RandSeed 1010 ► Done 
sample 

> {66. ,89. ,58. ,50. ,77. ,62. ,87. ,68. ,69. ,79. ,75. ,72. ,50. ,59. ,79. ,83. ,67. ,56., 47., 33. ,41. ,84., 65., 54. ,76. ,34. ,69. ,63. ,87. ,63. ,82., 6. ,62. ,87. ,61. ,80., 16. ,61., 82. ,6. ,44., 4 

samplemean=mean (sample) * 63.0938 

Number of calculated sample means: dim(means) * 1382 „ 

Compare the population mean and the mean of the sample means: 

mean(h4) ► 65.0027 mean(means) * 65.0682 

Compare the standard deviation of the population and the standard deviation of the sample means: 

stDevPop(h4) 


stDevPop(li4) ► 18.095 


stDevSampfmeaiis I 


2.21888 




2.26188 



• sample 

• means 

" 


=randsamp('h4, , n) 

=capture('samplemean,1 ) 



70 

63.5625 



77 

67.6563 



28 

63.6563 



77 

65.8281 



52 

61.3125 




60.8594 




63.8438 

i 

I 


68.1875 

I 



64.6563 




66.8125 




65.6875 




63.375 




66.75 


| A66 I 

-“m 


0.16- 


0 . 12 - 


0.08- 


0.04- 


0 . 00 - 


normPdf(x, 65. 0682, 2. 21888) 


25 


30 


35 


40 


45 


50 



75 


80 


85 


Sample size n = 64, 1382 samples. Sample is ... Illustration of the Central Limit Theorem. 


The reader is invited to experiment with samples taken from other (discrete or continuous) popula- 
tions. This shall lead to the conjecture that the probability model of the sample mean population - on 
condition that the sample size is sufficient large - can be approximated by a bell-shaped model - the 
normal distribution. 


The bell-shaped probability model of a normal distributed population is completely defined by its 
population mean // and its standard deviation a. The influence of these parameters on the graph of the 
density function will be investigated by introducing sliders. 

How to achieve the screen given below? This is the answer: 


• Split the page horizontally into a Graphs-application and a Notes-Application. 


• Activate the Graphs and introduce four sliders via the Graphs & Geometry Toolbox: 
1: Actions > B:Insert Slider. The Greek letters can be found among the Utilities in the 
Documents Toolbox. 


(Take care that you take the right p. 

The p-character shown in the first row is a 
system variable which cannot be used.) 

Insert sliders for two means fi u ju 2 and two 
standard deviation o x and cr 2 . 
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• Define the functions fl(x) and f2(x) via the entry line which appears after clicking on the 
» symbol bottom left (version 3.1). In version 3.2 make a right mouse click on the graph 
page. 2:Hide/Show > 5: Show Entry Line. 

• Define two points on the x-axis using Points & Lines > Point On in Graphs & Geometry 
Toolbox > 8:Geometry. Then press ESC for leaving this menu. 

• Show the coordinates of these points via 1: Actions > 8: Coordinates and Equations. 

• Switch to the Notes and define the variables a and b as shown below. 

• Return to the Graphs-application, move the cursor to the x-coordinate of the first point 
(the word “text” appears). Make a right mouse click on it and choose in the appearing 
pop-up menu 6:Variables > 2:Link to > a. Now this point is linked with the variable a. 
Repeat this procedure for the second point linking its x-coordinate with variable b. 


• Calculate the area under the graph of f2(x) between a and b using the Analyse Graph 
menu > 7:Integral. Click on the graph and then fix the two points as boundaries. 



Hint: In the Calculator Application and in the Spreadsheet Application you have access to the Maths 
Operators. Go to Statistics > Distributions and there you will find a couple of probability distribu- 
tions offered with the normal distribution among them. 

If X ~ N(p,o), then 

—> normPdf(x,|j 5 a) = f(x) with f being the normal density function 
-» normCdf(a,b,p,a) = P(tf<x<&) and 


-> invnorm(a 5 |j 5 a) = x with a = P(X < x) 


p 10 
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(5) Variability of samples in box plots 


It is worth while to observe the variability of the samples by generating box plots for samples follow- 
ing in a row. 


When taking small samples then the box plots should vary very much and should also differ from the 
population box plot significantly. Increasing the sample size the box plots of samples following each 
other are differing less and less. The random generated box plots shall also become more and more 
similar to the population box plot. 

Take again h4 (367 life guards data) as population and take random samples of sample size n. 


sample size n:=50 * 50 
sample tvew:-raiidSainp(h4,nJ 

* { 83, 72,70,75,6,42,84, 7 3,75,79,83, 80, 79,70,29, 8 2, 52, 14,90, 57, 50, 35,5 3, 90, 6 3,94,99, 58, 88,70,75, 81 , 75,78, 59, 84, 68, 5C 

I 

Simulate various samplejnews and observe their variability. 

Compare with the fixed box plot of the population h4. 


■'T 
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Random generated box plots with version 3.1 (above) and version 3.2 (below). 



Next DNL: Part 4: Discovering Probability Distributions 






D-N-L#87 


Surfaces from the Newspaper (7) 


P 11 


It was in DNL#71/72 when I presented some surfaces whose equations I found in 
the FAZ (Frankfurter Allgemeine Zeitung) some time ago. One empty page in this 
DNL is reason enough to have another three solids with DPGraph (left) and DERIVE 
(right), Josef 




r 2 2 2 3 ” 

> 



(x+y+z-1) 2 2 3 


VECTOR 

ContourDots_XY 

= 0 + y )-z , 1, 2, 2, 2, 2, 0.1, 0.1 

, 1 -1, 3, 0.1 



L io , 

j 



2 3 5 

VECTOR (Contour Pts_XY(x + y + z = 1, 


- 2 , 2 , - 2 , 2 , 0 . 1 , 0 . 1 ), 1 , - 2 , 2 , 0 . 1 ) 
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Overcoming Branch & Bound by Simulation 

Josef Bohm 

When I was a teacher at the St.Polten College for Business Administration (Handelsakademie) I had to 
teach several years “ Planungsmathematik ” which is a possible German word for Operations Research. 
This was in the early 70ties. The curriculum covered topics as Linear Programming, Critical Path 
Methods, and heuristic methods as Branch & Bound. Using Branch & Bound we solved Travelling 
Salesman problems, Knapsack problems and Assignment problems. Very small problems can be 
solved by complete enumeration, i.e. trying all possible combinations and checking them for validity 
and for leading to a maximum or minimum of the given goal function. 

All these problems are NP hard (they have no polynomial algorithm for solving). Bigger problems 
become very time consuming. 

Branch and Bound (B&B) is by far the most widely used tool for solving large 
scale /VP-hard combinatorial optimization problems. B&B is, however, an al- 
gorithm paradigm, which has to be filled out for each specific problem type, 
and numerous choices for each of the components exist. Even then, princi- 
ples for the design of efficient B&B algorithms have emerged over the 
years 1 * 1 . 

Guido Herweyers explains simulation performed with TI-Nspire. So I believe that it might be a good 
idea to demonstrate applying simulation for solving the problems mentioned above by simulation. 


1 The Travelling Salesman Problem ( Rundreiseproblem ) 

Given a set of locations, and known distances between each pair of loca- 
tions, the Travelling Salesman Problem is the problem of finding a tour that 
visits each location exactly once and that minimises the total distance trav- 
elled. 

Let’s start with 5 towns. The distances are given in the table below. The matrix needs not to be sym- 
metric. 


A 

A 

B 5 
C 10 
D 6 
E 7 


B C 
5 10 

8 

8 

8 8 
5 5 


D E 

6 7 

8 5 

8 5 

7 

7 


It needs 5 miles from A to B, 10 miles from A to C, etc. We will find the shortest trip starting in A 
visiting all other places and then returning to A. 

Route A-B-C-D-E-A would be 35 miles. 

[1] http://www.diku.dk/OLD/undervisning/2003e/datV-optimer/JensClausenNoter.pdf 
http://www.informatik.uni-trier.de/~ley /db/indices/a-tree/c/Clausen: Jens.html 
http://citeseerx.ist.psu.edu/viewdoc/summary7doU10. 1 . 1 .5.7475 
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But for now I will proceed on my travel route. I append the starting point because I have to finish the 
roundtrip and calculate the mileage which is again 35. Is this the shortest roundtrip? 


#10: 

[A, 

, c, 

D. E, 

B] 


#11: 

[A, 

, c, 

D, E, 

B, A] 


#12: 

di; 

stances 

1,3 

+ distances + dist; 

3,4 

inces + distances + distances = 35 

4,5 5,2 2,1 


It would be very boring to do this one route after the other. The solution is collecting all steps in one 
program: 

travel (dist, Iocs, numb, p, cost, way, dummy) := 

Prog 

[dummy := RANDOM(O), simul := []] 

Loop 

p := (SORT (VECTOR ([ RANDOM (1) , k], k, DIM(locs))))u2 
p := APPEND(p , [prl ]) 

#13: way := VECTOR ( locsj, k , k, p) 

cost := I(dist|(pr j)t(pi(j + 1)), j, 1, DIM(locs)) 
way := APPEND(way, [cost]) 
simul := APPEND(simul , [way]) 
numb 1 
If numb = 0 
RETURN simul 

I didn’t include simul in the parameter list because it is a global variable which is necessary for my 
intention./?, cost , way, and dummy are local variables. 

This is my first run of the program performing 5 simulations. Of course, It can happen that tours are 
appearing twice or even more often. 

#14: testrun := travel (distances , towns, 5) 

~ A B E C D A 29 ~ 

D B A C E D 35 

#15: testrun := B C E A D B 34 

E B C D A E 34 

_ A C D E B A 35 _ 

The first route is the shortest one with length 29. 

We can assume that the lengths L of the routes follow a normal distribution [1] . First we calculate mean 
and standard deviation of the sample and then the probability p{L < 29). 

#16: [AVERAGE (test run.u 7) , STDEV(testrunu?) ] = [33.4, 2.50998] 

f 29 - 33.4 ^ 

#17: NORMAL = 0.0398008 

L 2.50998 J 

There is a probability of only 4% to find a better solution, i.e. a shorter route combining all towns. 

A sample of size 5 is too small. I increase the number of simulations 1000 random routes. 
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#is 

#19 
#20 
#21 

#22 

When evaluating expression #19, the resulting matrix (1000 rows) is stored as global variable si mul . 

I need si mul for determining mean and standard deviation of the route lengths. If we would not take 
si mul in #21 but runl instead then a new simulation would be performed. If you trydim(simul) = 
you will get immediately 1000 but when simplifying di m(runl)= it will need some time to first per- 
form internally the next 1000 simulations and then show the number of rows = 1000. 

The probability to find a better solution has decreased to 2.5% which is pretty good. 

DIM(5ELECT (v = 29, v, simul)) = 87 
7 

~ A D C E B A 29 

(SELECT (v = 29, v, simul)) = CEBADC29 

7 

[1, 2, 3] L C E B A D C 29 

There are 87 routes with minimal length 87 which need not be all the same. I’d like to see the first 
three of them and am finding out that they all are the same. Maybe that there is at least one another 
one among the remaining 84 solutions. But we cannot be absolutely sure that 29 is the shortest possi- 
ble route. It remains an approximation. 

Selecting the minimum together with returning the optimal route could also be included into the pro- 
gram travel . I did without in the first application but I will do this for the other problems. 

You will see that my programming style will improve stepwise. This could also be done with students. 
I remember that at similar occasions the students had ideas what should - and often really could - be 
improved 

Let’s have a problem with more locations. This is a problem I used as teacher demonstrating the 
Branch & Bound Method. As I said above this was in the early 70ties. I will show my original “prepa- 
ration” provided for the classroom, doing all calculations by hand. The distance matrix - which is now 
defined as a cost matrix - is not symmetric and there are 7 locations. 


#23: 


#24: 


runl := travel (distances , towns, 1000) 
minimum(run) := MIN(run 4 , 4 ,DIM(runj,l)) 
minimum (runl) = 29 

[AVERAGE (simul j, j, 7) , STDEV(simulu 7) ] = [34.376, 2.73503] 
^ 29 - 34.3S 


NORMAL 


= 0.0254208 
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The last step performed on the next page shows that the branch below cannot provide a better solution 
than 107. 


The optimal route is A-B-C-D-G-F-E-A with minimal cost 107. 
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The DERIVE procedure: 


#27: 


#2S: 

#29: 

#30: 

#31: 

#32: 

#33: 

#34: 



CC' 

20 

33 

23 

12 

IS 

16 " 


19 

00 

24 

31 

20 

14 

20 


35 

25 

00 

26 

23 

22 

IS 

costm := 

24 

28 

28 

00 

16 

15 

11 


13 

23 

22 

17 

00 

6 

5 


20 

13 

22 

IS 

7 

00 

6 


_ 19 

19 

17 

12 

6 

6 

00 

offs := [A, B, 

c. 

D, E 

, F, 

G] 




run2 := travel (costm , offs, 1000) 
minimum(run2) = 107 

sel(val) := 5ELECT(v = val , v, simul) 

DIM(simulrl) 


sel (107) = [[A, B, C, D, G, F, E, A, 107]] 

[AVERAGE (simul 4 . - 1 - 9 ) , STDEV(simulu9) ] = [128.66, S. 33107] 


r 107 - 12S.66 ^ 


NORMAL 


S.53 


0.00535413 


The selection of all routes with a given value (107) is done by a function sel (val ). Comparing with 
my solution found in the last century we can see the correspondence. 
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The probability of 0.5% finding a shorter path is sufficiently small to accept 107 miles at least as a 
good approximation - even without knowing the exact result - as we do here. 

Finally let me show a more extended travel route as given in [2]: 

This is the map with the distances between villages A to M. We need to find the shortest path connect- 
ing all places starting and ending in town A. 

The first problem is to find out the 
shortest connections in pairs be- 
tween all villages - considering all 
nodes. We will assume that this 
problem has been solved for using 
a special algorithm (Dantzig). 

The distance matrix trm is given. 


6.5 12 14.2 17.2 12.2 11.2 17.8 24.1 20.2 15.5 9.2 

6.5 _ 5.5 7.8 10.7 5.7 6 12.2 18.2 18.6 12.8 12.7 

12 5.5 _ 6 14.1 11.2 11.5 17.7 20.2 24.1 19.2 19.2 

14.2 7.8 6 _ 8.1 6.5 12.8 12 15 . 6 20 . 8 21 . 6 21 . 5 

17.2 10.7 14.1 8.1 _ 5 16.7 11.5 7.5 19.2 22.5 24.4 

12.2 5.7 11.2 6.5 5 _ 11.7 6.5 12.5 14.2 18.5 19.4 

11.2 6 11.5 12.8 16.7 11.7 _ 6.5 12.8 12.6 7.8 9.5 

17.8 12.2 17.7 12 11.5 6.5 6.5 _ 6.2 7.8 12 16 

24.1 18.2 20.2 15.6 7.5 12.5 12.8 6.2 _ 12 18.2 22.2 

20.2 18.6 24.1 20.8 19.2 14.2 12.6 7.8 12 _ 6.5 18.5 

15.5 12.8 19.2 21.6 22.5 18.5 7.8 12 18.2 6.5 _ 12 

9.2 12.7 19.2 21.5 24.4 19.4 9.5 16 22.2 18.5 12 

[A, B, C, D, E, F, G, H, J , K, L, M] 
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1 1 villages give 1 814 400 possible routes. Which of them is the shortest one? I run three time 10000 
simulations: 


r un3 := travel (trm, vills, 10000)1 


#39: minimum(run3) = 111.5 

#40: 1500 sec 

#41: minimum(run3) = 111 

#42: minimum(run3) = 104.5 

#43: se] (104 . 5) = [ [C , B, E, J, H, L, K, G, M, A, F, D, C, 104.5]] 

#44: [AVERAGE(simu]ul4), STDEV(simulul4) ] = [165.684, 15 . 7818] 


#45: 


r 104.5 - 165. 68 ^ 


NORMAL 


15.78 


0.0000528632 


The best result is 104.5 miles. Calculation gives a very small chance to find a better one. 
Starting in A the route is:A-F-D-C-B-E-J-H-L-K-G-M-A. 

Preparing the DNL I tried once more hoping to find a better approximation: 


#47: minimum(run3) = 100.1 

#48: [ AVERAGE (s-imuiu^), STDEV(simu]ul4)] = [165.433, 15.6341] 

f 100.1 - 165.43 ^ 


#49: NORMAL, 

15.63 

#50: se] (100.1) = [ 


0.0000145894 


[E, F, H, J, K, L, G, B, A, M, C, D, E, 100.1] 


] 


The approximation given in [2] is much better. After some pages following a tiresome algorithm the 
approximation leads to a total distance of 84.1 when using the route A-B-C-D-F-E-J-H-K-L-G-M-A. 

The first figure shows the graph connected with the route length 104.5. 
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The next one gives the route with total length 84. 1 . 


Even be inspection we recognize that this is a much better solution. So we would need more simula- 
tions to obtain a shorter path? Maybe that we would be lucky with the next single simulation. The 
reader is invited to follow the route with length 100.1. 

The Travelling Salesman problem tackled with the TI-92 / Voyage 200: 

travel ( di s ,1 oc , numb) 

Prgm 

Local i , p ,c0 , cost , way , 1 i 1 , route 
ClrlO 
{ }->si ms 

X( di s [j , j + 1 ] , j ,1,dim(dis)[1]-1 )-»c0 
c0+di s [di m (di s ) [1 ] ,1 ]-»c0 
For i , 1 , numb 

seq( k,k,1 , di m(l oc) )-*p: seq ( rand ( ) , k, 1 , di m ( 1 oc) )-»l i 1 
SortA 1 i 1 , p :augment ( p , {p [1 ] } )-»p 
I(di s[p[j ] ,p[j + 1] ] ,j ,1 ,dim(p)-1 )-»cost 
augment ( si ms , {cost } )-»si ms 
If cost<c0 Then 
cost-»c0 : p-»way 
Endlf 
EndFor 

seq( 1 oc [way [i ] ] ,i , 1 , di m (way ) )-»way 

Di sp "Approximation for the cheapest route = " 

Di sp way 
Di sp " " 

Di sp "Cost = " &stri ng ( c0 ) 

EndPrgm 
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As you can see it works pretty well. 



raks. 


' FE Y f-.l* 

Prgn 1 0 |C L ::■<:! ■s 

4 I 

L7 5 5 

7 

L7 5 5 

7 

« C'FT "B" 

"C" "D" 

"E"} * towns 



£"fl" M E 

II II 011 II 011 

"E" y 

■ traueh! dist, towns, 50) 



■ nean^sims) 


34.720000 

■ stdDeu^sins) 

2.571746 

■ tistat. normpdf (29 , 34. 72 , 2. 57174601215fr 



.013076 


MAIN 

RAD AUTU 

DE H/ b 



The route produced by this simulation (50 tries) gives the optimal solution. 


|V : i!!!!!:]i” : vii'-Isr. 

.t'.'r.WS. 

'■fprghiof:; l sJr-T ] 

Approximation for the 

cheapest route = 

t"fl" "D" 
Cost = 29 

ii 0 ii 110-n 

M0M 

MAIN 

RAD AUTD 

DE H/ b 


I tried the second problem - with 7 locations: 


fit Y fjt y fht y fe y fat y ) 
|fllgebra|Calc|0ther|PrgmI0|dean Up| 



35 

25 

w 

26 

23 

22 

13 

■ dist2 


24 

23 

23 

« 

16 

15 

11 



13 

13 

22 

17 

w 

6 

5 



20 

13 

22 

13 

7 

w 

6 



.19 

19 

17 

12 

6 

6 

w 

■ Iocs 









{"FT 

ii 0ii 

"C" 

11 D 

II 

II 0 II 

"F 1 

II 

M 6" y 

t ravel Cdist 2, Iocs 

_300> 




MAIN 

RAD AUTD 


DE 1E/30 





VMu-iiua 1 

Approximation for the 

cheapest route = 

jT II p II II 0 II II 0 II II 0 II 

Cost = 103 

II 011 II p II II 011 II 0! 

MAIN RAD AUTD 

DE 13/30 


Not so bad! ! 


[1] Autorenkollektiv, Mathematische Standardmodelle der Operationsfors chung, Verlag Die 
Wirtschaft Berlin 1972. 

[2] Autorenkollektiv, Mathematische Standardmodelle der Operationsforschung, Aufgabensamm- 
lung, Verlag Die Wirtschaft Berlin 1972. 

http://en.wikipedia.org/wiki/Travelling salesman problem 
http://de.wikipedia.org/wiki/Problem des Handlungsreisenden 
http://www.tsp.gatech.edu/ 

http://mathworld.wolfram.com/TravelingSalesmanProblem.html 

http://www.lancs.ac.uk/staff/letchfoa/talks/TSP.pdf 

http://elib.zib.de/pub/UserHome/Groetschel/Spektrum/index2.html 


P 22 


Josef Bohm: Overcoming Branch & Bound 


D-N-L#87 


2 The Assignment Problem {Zuordnungsproblem) 

There are a number of agents and a number of tasks. Any agent can be 
assigned to perform any task, incurring some cost that may vary depending 
on the agent-task assignment. It is required to perform all tasks by assigning 
exactly one agent to each task in such a way that the total cost of the 
assignment is minimized. 

assign(ass, n, k, d, rows, cols, cn, rn, val , mx , mn, mxv, mnv, dummy) 
Prog 

[dummy := RANDOM(O), k := 1, d := DIM(ass), simul := []] 

[cn := REST(assil), rn := REST (ass 1 4,1) ] 

ass := VECTOR (VECTOR (ass|i_|j_, j_, 2, d) , i_, 2, d) 

d := DIM(ass) 

Loop 

If k > n exit 

rows := VECTOR (k_, k_, d) 

#1: cols := (SORT (VECTOR ([ RANDOM (1), k_] , k_, d)))||2 

val := I(ass| i| (colsi i) , i, 1, d) 
simul := APPEND(simul , [APPEND(cols , [val])]) 
k : + 1 

[mn := MIN(simul.n (d + 1)), mx := MAX(simulii (d + 1))] 
mnv := (SELECT (vr (d + 1) = mn, v, simul))il 

mxv := (SELECT (vr (d + 1) = mx , v, simul)) 4 ,l 

mn := [ 1 1 M I N : 11 , [rn, VECTGR(cn.j,k_, k_, mnv|[l d])], mn] 

mx := ["MAX: 11 , [rn, VECT0R(cn4,k_, k_, mxv|[l d])], mx] 

[mn; mx ; "all results in simul 11 ] 


The cost matrix - 6 projects are to be assigned to 6 agents is given by: 



PI 

P2 

P3 

P4 

P5 

i — 

kD 

Cl. 

A1 

10 

15 

S 

24 

12 

5 

A2 

4 

20 

6 

10 

S 

5 

A3 

20 

24 

10 

15 

12 

12 

A4 

15 

20 

12 

15 

10 

14 

A5 

30 

25 

20 

20 

25 

15 

_ A6 

12 

15 

15 

12 

15 

10 _ 


We can also interpret the matrix as a profit matrix because the agents perform their tasks more or less 
efficient. Then the maximum assignment is required. My program assi gn (mat ri x , numbe r of 
si mul ati ons) is looking for minimum and maximum as well. 

Additionally I have included all steps leading to the final answers (fixing min and max, selecting the 
respective assignments and providing a nicely formatted output of the results). 





" A1 A2 A3 A4 A5 A6 " 




MIN: , 

_ P5 P3 P6 PI P2 P4 _ 

, S2 

#3: assign(assl, 10) = 


MAX: , 

" A1 A2 A3 A4 A5 A6 " 

_ P6 P5 P2 P4 PI P3 _ 

, 97 


all results in simul 
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Having performed 10 simulations which will very likely not give good approximations, I will try 2000 
randomly chosen assignments. This needs only 9.7 seconds: 


#4: assign(assl, 200 


A1 

A2 

A3 

A4 

A5 

A6 

P6 

PI 

P3 

P5 

P4 

P2 

A1 

A2 

A3 

A4 

A5 

A6 

P4 

P2 

PI 

P6 

P5 

P3 

all 

resi 

nits 

in 

simu 

1 


#5 : needs 9 . 7 sec 

This is much better now. 

The probabilities to find better solutions are very small. 

#6: [ AVERAGE (simulul^), STDEVCsimululS)] 

#7: [86 . 7745 , 9.074074011] 

f 64 - 86.77 A 

#8: NORMAL = 0.006028458473 

L 9.07 J 

( 118 - 86.77 A 

#9: 1 - NORMAL = 0.000287416377 

L 9.07 J 


#9: 1 - NORMAL 


= 0.006028458473 


0.0002874163774 


Here is a second example from the 70ties-textbook. Seven jobs shall be assigned to seven machines. 
The cost/profit-matrix is given by: 

J1 J2 J3 J4 J5 J6 J7 
Ml 8 5 12 11 3 9 8 

M2 9 4 11 10 7 10 8 

M3 7 5 11 8 6 12 11 

ass2 := 

M4 8 8 9 7 2 10 7 

M5 5 7 11 7 8 8 11 

M6 11 2 7 5 6 8 7 

_ M7 7 5 2 10 8 11 10 


Looking back to my 40 years old papers I see that there are more optimal assignments possible. This 
reminiscence gives the opportunity for demonstrating the Branch & Bound method once more. 
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We see that the manual procedure gives four solutions for the minimum problem - and there are two 
solutions for the maximum problem. I changed my program as si gn in order to return not only one 
solution - but I cannot be sure that we will find all possible solutions. 
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#16: assign_more(ass2 , 10000) = 



" Ml 

M2 

M3 

M4 

M5 

M6 

M7 " 



6 

7 

2 

5 

1 

4 

3 


MIN: 

6 

7 

4 

5 

1 

2 

3 

36 


5 

2 

1 

7 

6 

4 

3 



_ 7 

2 

1 

5 

6 

4 

3 _ 



" Ml 

M2 

M3 

M4 

M5 

M6 

M7 " 



4 

3 

6 

2 

7 

1 

5 



3 

4 

6 

2 

7 

1 

5 


MAX: 

4 

3 

6 

2 

7 

1 

5 

72 


4 

3 

6 

2 

7 

1 

5 



_ 3 

4 

6 

2 

7 

1 

5 _ 



all results "in simul 


#17: 228 sec 


This took me 228 seconds. Comparing with my papers from the past, I can say that all optimal as- 
signments for minimum and maximum as well have been found by simulation. Look at the probabili- 
ties! 


[AVERAGE (simul 14 , 8 ) , STDEV(simuluS) 1 = [54.5964, 


NORMAL 


36 - 54.6 


5.66 


, 1 - NORMAL 


72 - 54.6 


5.66 


5.657005145] 

[0.0005077100026, 


0.001055321911] 


http://faculty.ksu.edu.sa/ikhan/Documents/ORyopt-lp5.pdf 
http://cse.spsu.edu/bmorrison/cs4413 ppt/branch&bound.ppt 



Josef, Karsten Schmidt and Carl Leinbach (Tartu 2012) 
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3 The Knapsack Problem ( Rucksackproblem ) 

Given a set of items, each with a weight and a value, determine the number 
of each item to include in a collection so that the total weight is less than or 
equal to a given limit and the total value is as large as possible. It derives its 
name from the problem faced by someone who is constrained by a fixed-size 
knapsack and must fill it with the most valuable items. 


Example: 

Let’s assume there is a burglar. He finds 6 items worthwhile to take them. They differ in weight and 
value. He is able to carry not more than 65 pounds. What is his best choice? 


Item / 

A1 

A2 

A3 

A4 

A5 

A6 

Weight w ■, 

30 

20 

8 

25 

12 

15 

Value V\ 

25 

18 

10 

20 

15 

8 


There is a recursive algorithm. The respective recursion formula is given by: 


m g ) 


\f(k-\,g) if w k >g 

jmax(v t + f(k-l,g- w k ), f(k - 1, g)) if w k < g and k> 0 


The last values are £ = (which is 6 in our case, and g = M= 65). 
The result of the procedure should be/(6,65) = 63. 


#1: simul := 

knapsack(supply , restr, n, i, z, art, amt, vals, m, r, ind, sel, sg, sv, mx , dummy) 
Prog 

[dummy := RANDOM(O), simul := []] 

i := 1 

[art := supplyj,l, amt := supplyr2, vals := supplyrB] 

[m := DIM(art) , mx := 1] 

Loop 

ind := VECT0R(k, k, m) 

If i > n ex it 
r := RANDOM (m - 1) + 2 
z := RANDOM (r) + 1 
sel := [z] 

#2: Loop 

If DIM(sel) = r exit 
ind := DELETE(ind, z) 

z := IF(DIM(ind) = 1, 1, RANDOM (DIM (ind)) + 1) 
sel := APPEND(sel, [ind.|.z]) 

[sg := I(VECTOR(amtrk , k, sel)), sv := I(VECT0R(vals4.k , k, sel))] 

If sg < restr a sv > mx 
Prog 

simul := APPEND(simul , [ [SORT (VECTOR (artr k , k, sel)), sg, sv]]) 
mx := sv 

i : + 1 

simul := SELECT(v|B = mx , v, simul) 

APPEND([ [’’Goods 11 , 'Total Weight 11 , 'Total Value 11 ]], simul) 
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See three series of simulations followed by the original binary B & B from my OR-teaching. 
#3: example := 


" A1 

A2 

A3 

A4 

A5 

A6 

30 

20 

S 

25 

12 

15 

_ 25 

IS 

10 

20 

15 

S 


#4: knapsack(example , 65, 10) = 


#5: knapsack(example , 65, 10) = 


Goods Total Weight Total Value 
[A1 , A3, A4] 63 55 

Goods Total Weight Total Value 


#6: knapsack(example , 65, 100) = 


[A2 , A3, AS, A6] 
Goods 

[A2 , A3, A4 , A5 ] 
[A2 , A3, A4 , A5 ] 
L [A2 , A3, A4 , A5 ] 


55 51 

Total Weight Total Value 


65 

65 

65 


63 

63 

63 
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This is a so called (0,l)-problem: Take it (= 1) or do not (= 0). So we have 2 6 possible choices from 
(0,0, 0,0, 0,0) - leave all items - to (1,1, 1,1, 1,1)- take them all. It is left to the reader to write an enu- 
merative procedure to check all possibilities (= end of all branches) for validity (total weight not more 
than 65 pounds) and for maximal total value. 


I transferred the Travelling Salesman problem to the TI-92 / Voyage 200. Now I will tackle the Knap- 
sack problem with TI-NspireCAS. 


The problem is from my old textbook again. It is a really “Knapsack = Rucksack” problem: 


A hiker packs his rucksack. He can to choose between eight victuals. 
Beef Pork Lard Butter Cheese Bread Chocolade Apples 

2 2 1 1 2 4 0.5 2 

_ 2600 S000 9200 S000 S200 10000 2300 1200 

The second row gives the weight (kg), the third one shows the calorie 
content. The hiker doesn’t want to carry more than 10 kg. 

Find the selection with maximum calorie content under the condition 
that the eatables are indivisible. 



This is the TI-Nspire program: 


knapsack 

Define knap s ack(sw/^ fy, restr,n)= 

Prgm 

Local i^simu^artfamtiValsJndfm^selfSw^sartfmx^ols 
szffzzz/:=[" Goods:" "Total Weight" "Total Value"] 
sols\=simul 

a rl : =mat ► 1 ist {supp />( 1 ] ) : a m t : =mat ► 1 ist {supp />( 2 ] ) : va Is : =mat ► 1 ist {supp />[ 3 ] )| 
zrcx: = l:ffz:=dim(tfz’f) 


For i } l } ti 

i n d : =seq ( k, k, 1 1 m ) : r : =randlnt ( 2j m ) 
se 1 : =rand Samp ( i n d t r, 1 ) : SortA sel 

dim(se/) dim(se/) 




j> (va/s[se/[?]]) 


J = * 1 J= l 

sa rl : =str in g (seq ( a [se /[/ ] L 1 , d im(se /) ) ) 

If sw<restr and s\>>mx Then 

si m it I : =cg lAu gmentfsz m u ! } [sa rl sw sv]):/77.x: =sv 

Endlf 


EndFor 

For ^2^dim(s7mu/)[l] 

If simul[i } 3\=mx Then 

so Is : =co lAu gment (so Is, si mu /[ i ] ) 

Endlf 


EndFor 
Disp sols 
EndPrgm 
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"{"A1","A3","A5'7'A6"}" 

"{"A1" J M A2";A5"}" 


knapsac^example, 65, 1 oo) 




"Goods:" 


"Total Weight" 

"Total Value" 


’> 

to 

"A3" J I! A4" J "A5 I 

}" 

65 

63 


"{"A2", 

"A3" J "A4" J "A5 I 


65 

63 


"{"A2", 

"A3" J II A4" J "A5 I 


65 

63 


"Beef" "Pork" "Lard" "Butter" M Cheese" "Bread" "Chaco lade 11 "Apples" 
2 2 1 1 2 4 0.5 2 

2600 8000 9200 8000 8200 10000 2300 1200 



"Beef" "Pork" 

"Lard 

11 "Butter" 

"Cheese" "Bread" 

"Chocolad' 

e" "Apples" 


2 2 

1 

1 

2 4 

0.5 

2 


_ 2600 8000 

9200 

8000 

8200 10000 

2300 

1200 

knapsack{hike } 10,100) 



"Goods:" 

"Total Weight" 1 

'Total Value" 

" {"Beef 1 

/’Pork", 

"Lard", "Butter", "Bread"}" 

10 

37800 


knapsack{hike t 10, loo) 


"Goods:" "Total Weight" "Total Value 

"{"Beef ', "Lard", "Butter", "Cheese", "Bread" }" 10 38000 


knapsack{hike ,10,500) 


"Goods:" "Total Weight" "Total Value 

"{ "Pork", "Lard", "Butter", "Cheese", "Bread"} 11 10 43400 

"{"Pork", "Lard", "Butter", "Cheese", "Bread"}" 10 43400 


And what is DERIVE recommending? 


knapsackfhike , 10, 100) = 


knapsackfhike , 10, 500) = 


Goods Total Weight Total Value 

[Beef, Bread, Butter, Lard, Pork] 10 37S00 

[Beef, Bread, Butter, Lard, Pork] 10 37800 

Goods Total Weight Total Value 

[Bread, Butter, Cheese, Lard, Pork] 10 43400 


Following the result of the first mn with only 100 simulations performed we would miss 5600 calo- 
ries. Running 500 simulations needs only 1 second and gives a much better result. 
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knapsack(supp]y , restr, n, i, z, art, amt, vals, m, r, ind, sel, sg, sv, mx , dummy) 
Prog 

[dummy := RANDOM(O), simul := []] 

i := 1 

[art := supply|l, amt := supplyj.2, vals := supplyi.3] 

[m := DIM(art) , mx := 0] 

Loop 

-ind := VECT0R(k, k, m) 

If i > n exit 
r := RANDOM (m - 1) + 2 
z := RANDOM (r) + 1 
sel := [z] 

#1 : Loop 

If DIM(sel) = r exit 
"ind := DELETE(-ind, z) 

z := IF(DIM(ind) = 1, 1, RANDOM (DIM (ind)) + 1) 
sel := APPEND(sel, [ind^z]) 

[sg := I(VECT0R(amtik , k, sel)), sv := I(VECT0R(vals|k , k, sel))] 

If sg < restr a sv > mx 
Prog 

simul := APPEND(simul , [ [S0RT(VECT0R(art|k , k, sel)), sg, sv]]) 
mx := sv 

i : + 1 

simul := 5ELECT(v|3 = mx , v, simul) 

APPEND([ [’’Goods 11 , ’’Total Weight”, ’’Total Value' 1 ]], simu]) 


I found in Muller-Merbach [1] an economic application: An investor wants to invest 1300 monetary 
units (say 1000 €). He has the choice of 7 indivisible investments with different returns of investment 
(also given in 1000 €). How should he put together his portfolio? 



" A 

B 

C 

D 

#6: stocks := 

500 

□ 

o 

LTl 

400 

300 


_ 35 

31 

26 

20 


#7: knapsack(stocks , 1300, 50) = 


#S: knapsack(stocks , 1300, 100) = 


#9: knapsack(stocks , 1300, 300) = 


(knapsack(stocks , 1300, 1000)) 

# 10 : 


E F G " 

200 100 100 
IS S 6 _ 

Goods Total Weight Total Value 


[A, B, D] 1300 S5 J 

Goods Total Weight Total Value 

[A, B, E, F] 1300 92 

Goods Total Weight Total Value 

[A, C, E, F, G] 1300 93 

[A, C, E, F, G] 1300 93 

[A, C, E, F, G] 1300 93 


= [[A, C, E, F, G], 1300, 93] 


Our recommendation is: buy stocks A, C, E, F and G and make a profit of 93 000 €! ! 





D 
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[1] Muller-Merbach, Operations Research , Vahlen 1971 

The Operations Research Problem Solver , REA, New York 1985 

http://www.diku.dk/users/pisinger/95-l.pdf 

http://en.wikipedia.org/wiki/Knapsack problem 

http://www.es.ele.tue.nl/education/5MC10/Solutions/knapsack.pdf 

http://www.diku.dk/users/pisinger/95-l.pdf 

http://www.nuigalway.ie/mat/algorithms/knapsack.html 
(a program!!) 

http://www.cs.ucsb.edu/~iacopo/BB/ 

(a Branch & Bound Solver) 

I’d like to add two comments: 

(1) I really like to undertake extended hikes. But I never had the idea to carry lard with me. It 
might be difficult enough to have butter in the rucksack climbing up a mountain on its sunny 
side in summer time .... And I never had beef and pork with me. And why should I not be able 
to divide butter, bread, cheese, . . .? Believe it or not this is the original text from the textbook. 

(2) I also like programming and exploring the advantages and disadvantages of various program- 
ming languages and various programmable computer algebra systems. I missed Nspire’s 
randSamp function. So what to do? Easy answer: make your own DERIVE randSamp- 
function and benefit from DERIVE’s feature to implement “real” random processes without 
needing a “randseed” using its RANDOM(O) command. 



4 randSamp & co 

Let’s start with a random permutation of n elements. I wantd to avoid simplifying RANDOM(O) at the 
beginning of each session so I included this command in the function in order to have every time an- 
other initial value for the random number generating process. 


#1: test := [A, B, C, D, E, F, G, H] 

randp(l, dummy) := 

Prog 

#2: dummy := RANDOM (0) 

VECTOR ( 1 j, k_, k_, (SORT (VECTOR ([RAN DOM (1) , k], k, DIM(l))))u2) 

#3: randp(test) = [E, A, H, G, C, F, D, B] 

This looked pretty good. I proceeded with a sequence of three random permutations of the set test in 
a row. 


#4: 


VECTGR(randp(test) , k, 3) 


"ABDFECGH" 

ABDFECGH 

.ABDFECGH. 


I was very disappointed that the same permutation is appearing three times! 
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Ok, I’ll come back to this problem later. We wanted to have random samples - with and without re- 
placement (of the elements). 


#6: 

randsamp(test , 

5) = [F, 

A, G, 

E, 

E] 


#7: 

randsamp(test , 

5) = [E, 

H, F, 

E, 

F] 


#&: 

randsamp(test , 

5,1) = 

l 1 

M 

O 

H, 

B, 

A] 

#9: 

randsamp(test , 

5,1) = 

i — i 

m 

T| 

G, 

B, 

C] 

j 

#10: 

VECTOR(randsamp(test , 5 

, i), 

k, 

3) = 

■ 


" B 

A 

E 

H 

C " 

B 

A 

E 

H 

C 

_ B 

A 

E 

H 

C _ 


The same problem as above occurred. One function call after the other worked properly, a sequence of 
random samples seemed not to be really randomly! ! randsamp(popul ati on , n) generates a sample 
of n elements with repetition by default). Including 1 as third parameter avoids repetition of elements 
in the sample. 

Do you know the reason for the repeated samples? I believe to know. The CPU works too fast. RAN- 
DOM^) uses the internal CPU time as initial value for the random number generating algorithm. And 
it seems to be that this time does not change between two or even more calls of the function. What to 
do is, giving the machine some time for “thinking” or doing some useless work: 

randsamp(l, n, s := 0, dummy, "i, j, samp) := 

Prog 

dummy := RANDOM (0) 

i := 1 

samp := [] 

Loop 

If i > n 

RETURN samp 

#11: k := IF(DIM(1) > 1, RANDOM (DIM(l)) + 1, 1) 

samp := APPEND(samp, [Irk]) 

If s = 1 

1 := DELETEd, k) 

j := 0 

Loop 

If j = 100 exit 
J :+ 1 

i :+ 1 

#12: randperm(l) := randsamptl, DIM ( 1 ) , 1) 

I let DERIVE count from 1 to 100 in the last loop. Counting only to 10 or 20 is not sufficient (I tried). 

I can use randsamp to produce a random permutation, too, because this is only a special case (taking a 
sample of all elements without replacement). 


#13: 

randsamp(test , 

5) = [F, 

c, 

A, 

C, 

C] 

#14: 

randsamp(test , 

5,D = 

[D, 

B , 

F, 

i 1 

< 

#15: 

randperm(test) 

= [C, H, 

G, 

B, 

F, 

i — i 
□ 

< 

LU 


I try again generating a sequence of samples and permutations: 
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" F C B F E H " 

#16: VECTOR (randsairip (test, 6), k, 3) = DHHGDE 

_ D C B E B G _ 

" H B G F E C " 

#17: VECTOR(randsamp(test , 6, 1), k, 3) = DFGCEB 

_ D A B E F H _ 

"AGBHECFD" 
HDBEGCAF 
#18: VECTOR(randperm(test) , k, 5) = GAECBHFD 

GHDBEFCA 

_DEHCFBAG_ 

Now everything is as expected. See other sequences of samples: 

#19: VECTOR(randsamp(test , k) , k, 5) 

#20: [[H], [B, A], [B. A, H], [D. B, F, D] , [A, F, B. F. H]] 

#21: VECTOR(randsamp(test , k, 1), k, 5) 

#22: [[A], [E, B], [E, B, D] . [G, C, A, D] , [D. H, E, C. G]] 

#23: VECTOR(randsamp(test , k) , k, S) 

Speaking about permutations, I wanted to generate pairs of permutations where no single element 
remains on its position (permutations without a fix point). 

diffperms(l , pi, p2, i) := 

Prog 

pi := randperm(l) 
p2 := randperm(pl) 

Loop 

i := 1 

#29: Loop 

If plj,i = p2j.. "i exit 
If i = DIM(pl) 

RETURN [pi, p2] 

i :+ 1 

p2 := randperm(pl) 

"HGEFBACD" 

#30: diffperms(test) = 

Lfabecgdh. 


#31: VECTOR(diffperms(test) , k, 1, 2) 



FAHCEDBG 

EBFGDCHA 



#33: d i f f pe rms (VECTOR (k , k, 1, 13)) 

' 6 14 2 13 7 10 12 1 11 13 S 3 4 3 9 

# 34 : 

_ 1 3 10 3 14 11 7 4 2 9 13 6 13 S 12 





p 34 


Josef Bohm: Overcoming Branch & Bound 


D-N-L#87 


I intended to present a fourth Branch & Bound problem which needs the pairs which can be formed by 
numbers 1 to n without replacement and without considering the order (i.e. the combinations without 
repetitions) 


pa-irs(l) := APPENDf VECTOR (VECTOR (T1 , 1 1, j , i + 1, DIM(l)), i, 1, DIM(l) - 1)) 
#35: L i jJ 


" 12 " 
1 3 


#36: pan rs(VECTOR(k , k, 1, 4)) 


1 4 

2 3 


2 4 


L 3 4 J 


The number of possible pairs is given by • So we should have 190 pairs formed of 20 elements. 

#37: DIM(pai rs(VECTOR(k , k, 1, 20))) = i0g| 

I save this file as randsamp . mth for possible use as utility file. 


5 The Localisation Problem 

We have the problem to distribute n offices (workstations, departments, ...) 
among n rooms (factories, districts, towns, ...). The distances between the 
localities are given and the intensity of necessary communications between 
the offices as well (number of travels, duration of communication, ...). 

How to distribute the offices among the rooms in order to minimize the sum of 
all communications (travel km, time needed, ...)? 

Example: 

The distances r ki between rooms Rl, R2, R3, R4, and R5 (in minutes walking distance) are given in 
form of a matrix (r_val s) and the average number of communications by between the offices A, B, C, 
D, and E is given in matrix b_val s. 


r u 

Rl 

R2 

R3 

R4 

R5 

b , 

A B 

C 

D 

E 

Rl 

- 

2 

7 

5 

5 

A 

- 15 

20 

12 

14 

R2 


- 

3 

2 

6 

B 

- 

10 

15 

12 

R3 



- 

4 

4 

C 


- 

8 

10 

R4 




- 

3 

D 



- 

6 

R5 





- 

E 




- 


For programming reasons I need one of the symmetric matrices in full form, and I chose the distance 
matrix r_val s. 


D 
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" 0 

2 

7 

5 

5 " 



" 0 

15 

10 

12 

14 " 















2 

0 

3 

2 

6 



0 

0 

10 

15 

12 








b_val s := 






, r_vals := 

7 

3 

0 

4 

4 



0 

0 

0 

S 

10 















5 

2 

4 

0 

3 



_ 0 

0 

0 

0 

6 _ 














_ 5 

6 

4 

3 

0 _ 



Read this matrices as follows: 

15 communications between A and B if A is located in room 1 and B in room 3 will need totally 
b\ 2 *r 13 = 15 • 7 = 105 minutes, but if A and B were located in rooms 5 and 3 then the time needed for 
walking between A and B would be b\rr 53 = 15 • 4 = 60 minutes, etc. 

For understanding my program let’s assume that the present distribution is given by: 


Offices 

A 

B 

C 

D 

E 

Rooms 

R3 

R4 

R2 

R1 

R5 


All necessary walks sum up as follows to 467 minutes. 


(A,B) • (R3,R4) 

+ 

(A,C)(R3,R2) 

+ 

b 12‘ r 34 

+ 

A ‘ r 32 = r 23 

+ 

15-4 

+ 

10-3 

+ 



623 • 

+ 



10-2 

+ 


(A,D)(R3,R1) 

+ 

(A,E)(R3,R5) 


II 

U> 

+ 

As ’ ^35 


12-7 

+ 

14-4 

230 

^24 ' r i4 

+ 

b 25 ' r 45 


15-5 

+ 

12-3 

131 

‘ r \2 

+ 

b 35 ' r 25 


8-2 

+ 

10-6 

76 







6-5 

30 


467 

Now you can recognize the purpose of the list of pairs (named co in the following program). 

I preload randsamp . mth as utility file to use randperm and pai rs as well. 

#1 : L0AD(D: \D0KUS\DNL\DNLS7\ randsamp . mth) 

locprob(bv, rv, numb, d, p, 1, co, l_opt, cost, opt, i) := 

Prog 

[d := DIM(rv) , p := VECTORfk , k, 1 , d) , opt := 100000 ] 

[i := 1, co := pairs(p), simul := []] 

Loop 

If i > numb exit 
1 := randperm(p) 

cost := I(bvr (cor jri)r (cor ■ rvr (lr (cor jri))r(lr (cor jr^)) i j, l, dim(co)) 
#2: WRITE([i, cost]) 

simul := APPEND(simu] , [cost]) 

If cost < opt 
Prog 

opt := cost 
]_opt := ] 

i : + 1 

DISPLAY( n al 1 results in simul 11 ) 

[l_opt , APPEMD( n Min = 11 , cost)] 


The core line for finding the sum as described above is given by cost :=.... 
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I run 100 simulations: 


#4: locprob(b_vals , r_vals, 100) 

all results in simul 

#5: [[4, 2, 5, 1, 3], Min = 424] 

#6: [AVERAGE (simul) , STDEV(simul) ] = [459.3, 14.0023357] 

#7: NORMAL (424, 459. 3, 14.0023357) = 

The recommended distribution is: 


Offices 

A 

B 

C 

D 

E 

Rooms 

R4 

R2 

R5 

R1 

R3 


giving a total sum of 424 walking minutes per period. The probability to find a better solution is only 
0.5%. 

My paper & pen solution from 40 years ago is also 424 - but it needed some paper work! ! 


The program includes two little tricks which might not be well known to all DERIVIANs: 

(1) TheWRITE([i ,cost]) command presents the progress of the program. The number of the 
simulation just performed together with its result is given in the bottom left corner of the 
status line. 


[289, 440] 


Of 

P 

V 

5 

E 

A 

B 

r 

A 

E 


This is a part of a screenshot taken during a 500 simulations’ 
run. We see that we have performed 289 simulations so far 
and the cost of this distribution is 440 minutes. 

WRITE is very useful for debugging programs. 


(2) At the end of the program I used the DISPLAY -command instead of including the text into the 
program output. This can be useful if you want to give additional information - and also for 
testing during developing a program. 


I will close my expedition in the world of Branch & Bound with a more extended localisation prob- 
lem: Matrix acti vi ti es shows the number of activities between departments A to H and matrix 
ti mes gives the transport times between sites R1 to R8. How to distribute the departments? 











" 0 

13 

26 

37 

15 

3 

12 

10 " 



" 0 

3 

4 

5 

6 

2 

10 

7 " 





















13 

0 

9 

25 

31 

46 

10 

3 



0 

0 

1 

3 

5 

4 

2 

9 





















26 

9 

0 

34 

43 

19 

17 

5 



0 

0 

0 

11 

3 

12 

1 

3 





















37 

25 

34 

0 

6 

49 

22 

11 


activities := 

0 

0 

0 

0 

6 

5 

10 

0 

, times := 




















15 

31 

43 

6 

0 

53 

29 

16 



0 

0 

0 

0 

0 

7 

2 

6 





















3 

46 

19 

49 

53 

0 

17 

9 



0 

0 

0 

0 

0 

0 

9 

1 





















12 

10 

17 

22 

29 

17 

0 

20 



_ 0 

0 

0 

0 

0 

0 

0 

5 _ 




















_ 10 

3 

5 

11 

16 

9 

20 

0 _ 
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Josef Bohm: Overcoming Branch & Bound 
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#9: locprob(activit"ies , times, 1000) 

all results in simul 


#10: [[4, 2 , 6, S, 7, 1, 5, 3], Min = 2611] 

#11: [AVERAGE(simul) , STDEV(simul) ] = [3286.18, 240.7044416] 

#12: NORMAL (2611, 3286.18, 240.7044416) = 0.002313697936 


My approximation is: 


Departments 

A 

B 

C 

D 

E 

F 

G 

H 

Localities 

R4 

R2 

R6 

R8 

R7 

R1 

R5 

R3 


with a minimal amount of time which is 2611 time units. (You may double check the result!) 

There is tiny probability of 0.2% to find a better way to distribute the departments, so we are pretty 
fine with the result. 

I had much pleasure doing all the programming and I am quite sure that there are some ways to 
improve the procedures. I would like to invite you to improve the programs and to provide other B & 
B methods for simulation. 



Many DUG members in front of the Kissing 
Students on the town square of Tartu 
(above). 


Part of the TIME 2012 social program: 
Learning and practising counting in Estonian 
(right). 
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Comment on “Circles in a Point” (DNL#86) 

Lieber Herr Bohm, 

vielen Dank fur die neue Ausgabe des DNL ! 

Besonders gelungen finde ich diesmal den Artikel ^Points in a Circ- 
le" von Roland Schroder. 

Bei der Anpassung der Losung fur die TI-Nspire-Sof tware haben Sie 
mangels VECTOR-Bef ehl eine andere Losung vorgeschlagen . 

Mit dem Nspire-Bef ehl constructMat lasst sich Herrn Schroders Vorge- 
hensweise jedoch weitgehend identisch ubertragen: 

Many thanks for the new DNL issue. In particular I liked Roland Schoder’s contribution on 
“Points in a Circle” very much. 

Adapting the solution for Tl-Nspire you proposed another approach because you missed the 
VECTOR command. But by applying the Nspire command constructMat it is possible to 
transfer Roland’s procedure nearly identical. 


constructMat when 


r+ 1,2- r+l)-*g(r) 


Per tig 


f(4) 

0 

0 

0 

0 

1 

0 

0 

0 

0 



0 

0 

1 

1 

1 

1 

1 

0 

0 



0 

1 

1 

1 

1 

1 

1 

1 

0 



0 

1 

1 

1 

1 

1 

1 

1 

0 



1 

1 

1 

1 

1 

1 

1 

1 

1 



0 

1 

1 

1 

1 

1 

1 

1 

0 



0 

1 

1 

1 

1 

1 

1 

1 

0 



0 

0 

1 

1 

1 

1 

1 

0 

0 



_0 

0 

0 

0 

1 

0 

0 

0 

0 _ 


count!f(g(4),l) 
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Herzliche GriiBe / Best regards 

Jurgen Wagner, DUG Member since 1992 


The Tribonacci Sequence with Tl-Nspire 

In DNL#86 I did not use recursion for producing this sequence with TI-NspireCAS because 
I thought that this would be impossible - and I complained in Tartu about the lack of recur- 
sion. But Michel Beaudin gave a very valuable advice how to perform a recursion: 



1 , 

x=l Done 


lrib(x):=' 

1 , 

x =2 



2 , 

x =3 



trib[x-\\+tril 

)[x- 2 )+trib\x-?>\ ,x >3 


tribi l) 


1 


trib[ 4 ) 


4 


seq(rn’Z?(/c;l,/r,l,lo] 

{1,1,2,4,7,13,24,44,81,149} 



Many thanks to Michel. It is my pleasure to revise my opinion. Recursion is possible, Josef. 


D 
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DERIVE 4 for DOS on HP 200LX 

Minh Nguyen rmailto:mn1 993@q raffiti.net! 

Hi, I am looking for derive 4 for dos to run on my HP 200LX. you know where i can get this 
now? Can’t find it anywhere. Much appreciated. Thanks 

Two problems provided by Fred Tydeman and Karsten Schmidt 

Problem 1: Non Recursive Definition 

Fred Tydeman wrote on 1 1 August 2012: 

I have a set of values and two recursive definitions: 

/(l) = 0 

/( 2) = 2*0 + l = l*l = l 
/(3) = 3*l-l = 2*l = 2*(l + 0) = 2 
/(4) = 4*2 + l = 3*3 = 3*(2 + l) = 9 
/(5) = 5*9-1 = 4*1 l = 4*(9 + 2) = 44 
/(6) = 6* 44 + 1 = 5* 53 = 5* (44 + 9) = 265 


The generalisation reads as follows: 

= + /( 1 ) = 0 
or 

m = (n - 1) * (f(n - 1) + f{n - 2)); /( 1) = 0, /( 2) = 1 

Is there some way to get a non-recursive definition? 

His comment on 12 August was: 

Thanks to all that responded. They used either Maple or MATHEMATICA to solve this. 

No one used DERIVE. 

In looking at the various solutions sent to me, I came up with this: 

TL ! 

/(«) = — , rounded to the nearest integer. 
e 

This should be a challenge for the DERIVIANs - and the “Nspirians” as well to tackle 

this problem! 
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Problem 2: Eigendecomposition 

. . . and for this I need the eigendecomposition (spectral decomposition = Spektralzerlegung) of a 
symmetric matrix: 

Let A a symmetric matrix. Then A = QLQi L is a diagonal matrix with the eigenvalues of A in the 
main diagonal and Q is the orthogonal matrix of the eigenvectors of A (see rule 4.1.1 1 in our Springer 
textbook [1] ). 

Do you know if anybody has done this with DERIVE before? This is, you know, very tiresome to 
write the eigenvalues into the main diagonal of a diagonal matrix. (The eigenvalues can appear multi- 
ple, i.e. the number of them which are returned by the EIGENVALUES function can be too small.) 
And it is much more laborious to calculate the respective eigenvectors using the EX- 
ACTEIGENVECTORS-function, scale them to length 1, and then write them into matrix Q. 

Best regards 
Karsten Schmidt 


I am no matrix algebra specialist at all. I took Karsten’s textbook and studied carefully rule 
4.1 .1 1 presenting the spectral decomposition of a 2x2 matrix together with the following ex- 
ample 9. I wanted to apply this decomposition on a 3x3 example. On the web I found (I don’t 
know which CAS was used): 


Example: Spectral Decomposition of a Correlation Matrix [2] 

C <- matrix (nrow=3, ncol=3, data=0.5) 

C[l,2] = C[2,l] <- 0.7 

C[2,3] = C[3,2] <- 0.9 

diag(C) <- 1 

print (C) 


gives 

[,1] 

[ / 2 ] 

[ / 3 ] 

[1, ] 

1 . 0 

0.7 

0.5 

[2, ] 

0.7 

1 . 0 

0.9 

[3, ] 

0.5 

0.9 

1 . 0 


Eigenvalue Decomposition: 

sz <- eigen (C) 

# matrix of eigenvectors 
Gamma <- sz$vectors 

# diagonal matrix of eigenvectors 
Lambda <- diag ( sz$values ) 

# test, if it results as C again 
test <- Gamma %*% Lambda %*% t (Gamma) 
print (test) 


gives 

[,1] 

[ / 2 ] 

[ / 3 ] 

[l, ] 

1 . 0 

0.7 

0.5 

[2, ] 

0.7 

1 . 0 

0.9 

[3, ] 

0.5 

0.9 

1 . 0 
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I tried simplifying Q Q' but I could not wait for the result, so I approximated the product and 
was pretty sure to receive the unity matrix (#1 1) - more or less approximated - as expected. 
Finally I performed the “test” from above using the approximated matrices only. 


#12: Q ■ diag(ev) -Q 1 


1 0.7 0.5 ' 

0.7 1 0.9 

0.5 0.9 1 . 


I had expected a result close to the given matrix A because of my numerical procedure, but - 
and this was my second, but positive surprise - the result was exact matrix A (in decimal 
form!). To be honest, if you increase the number of decimal places you will recognize that the 
result is rounded. This is the story of my first successful eigendecomposition of a symmetric 
matrix. 

[11 Schmidt, Trenkler, Einfuhrung in die Moderne Matrix-Algebra, Springer 2006 


I tried with TI-NspireCAS. I did not benefit from diag() only but also from the fact that 
eigVc(matrix) immediately returns the orthogonal matrix of eigenvectors. At the other hand 
Tl-Nspire does not provide any exact results. 


ev:=eigVc 

| 1 0.7 0.5 

0.7 1 0.9 

J.0.5 0.9 1 _ 


0.516947 0.818814 0.249618 

0.627235 -0.163883 -0.761393 
0.582531 -0.550169 0.598308_ 


i 

dev:= diag 

eigVl 

1 0.7 0.5| 

0.7 1 0.9 

,0,5 0.9 1 Jj 


2.41277 0 0 

0 0.523942 0 

0 0 0.063283_ 


ev dev ev T 

1. 0.7 0.5 

0.7 1. 0.9 

0.5 0.9 1. _ 



Karsten sent two additional DERIVE files: 
spectral2x2 . dfw (treating example 9 from [1] 

spectral 3x3F (decomposition of a 3x3-matrix with multiple eigenvalues) 


http : // www . stat . uni-muenchen . de/~chris/ rkurs/Rkurs6 . pdf 


http : / /www . ece . ucsb . edu/~roy/ class notes/ 210a/ECE210a lecture 15 small . pdf 


http : / / cseweb . ucsd . edu/~dasgupta/2 91/lec7 . pdf 


http : / / www . youtube . com/ watch?v=hmahoeta4KY 


